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ABSTRACT 



We report on the morphological and spectral evolution of 
SNR 1987 A from the monitoring observations with the 
Chandra/AClS. As of 2005, the X-ray-bright lobes are 
continuously brightening and expanding all around the 
ring. The softening of the overall X-ray spectrum also 
continues. The X-ray lightcurve is particularly remark- 
able: i.e., the recent soft X-ray flux increase rate is sig- 
nificantly deviating from the model which successfully 
fits the earlier data, indicating even faster flux increase 
rate since early 2004 (day ^6200). We also report results 
from high resolution spectral analysis with deep Chan- 
dra/LETG observations. The high resolution X-ray line 
emission features unambiguously reveal that the X-ray 
emission of SNR 1987 A is originating primarily from 
a "disk" along the inner ring rather than from a spheri- 
cal shell. We present the ionization structures, elemental 
abundances, and the shock velocities of the X-ray emit- 
ting plasma. 
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1. INTRODUCTION 



The continuous development of optically bright spots 
around the inner ring of supernova (SN) 1987A indicates 
that the blast wave shock is approaching the dense cir- 
cumstellar material (CSM) that was produced by the mas- 
sive progenitor's equatorial stellar winds. In such a case, 
the blast wave shock front should decelerate by interact- 
ing with the dense CSM, and may produce significant soft 
X-ray emission. We have been monitoring the evolution 
of SN 1987 A with the Chandra X-Ray Observatory since 
1999, and found that the soft X-ray flux from SN 1987A 
has indeed been rapidly increasing (Burrows et al. 2000; 
Park et al. 2002;2004). The significant interaction of the 
blast wave with the dense CSM therefore signals the birth 
of a supernova remnant (SNR) 1987 A. We successfully 



described the composite soft X-ray (0.5—2 keV) light 
curve, obtained with the ROSAT and the Chandra be- 
tween 1990 and 2003, using a simple model that assumes 
a constant-velocity shock interacting with an exponential 
ambient density profile (Park et al. 2004;2005a). We re- 
port here the latest development of the X-ray light curves 
of SNR 1987 A, in which we find a significant deviation 
of the soft X-ray flux from the previous model predictions 
since day ^6200. 

In addition to the regular monitoring observations, we 
have performed two deep Chandra gratings observations 
of SNR 1987 A. The early observations with the High 
Energy Transmission Gratings Spectrometer (HETG) in- 
dicated a high-velocity shock (v ^ 3400 km s~^) in 
which the X-ray emitting plasma is in electron-ion non- 
equilibrium (Michael et al. 2002). SNR 1987A was, 
however, faint and the detected photon statistics were 
limited. Recently, we performed follow-up deep obser- 
vations of SNR 1987 A with Low Energy Transmission 
Gratings Spectrometer (LETG) (Zhekov et al. 2005a). 
With the SNR being ^8 times brighter and the exposure 
^^3 times deeper than the early HETG observations, the 
new LETG data provide good photon statistics and allow 
us, for the first time, to measure the individual line pro- 
files and line ratios from SNR 1987 A. We present here 
some first results from these LETG data. 



2. OBSERVATIONS 



As of 2005 July, we have performed a total of eleven 
monitoring observations of SNR 1987 A with Advanced 
CCD Imaging Spectrometer (ACIS) aboard Chandra (Ta- 
ble 1). We also have performed two deep gratings ob- 
servations in 1999 October and 2004 August-September 
(Table 1). The data reduction of these observations have 
been described in literatures (Burrows et al. 2000; Park 
et al. 2005c; Michael et al. 2002; Zhekov et al. 2005a). 
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Table 1. Chandra Observations ofSNR 1987A. 



Observation ID 


Date (Age^) 


Instrument 


Exposure (ks) 


Counts 


124+1387^ 


1999-10-6(4609) 


ACIS-S + HETG 


116.1 


690^ 


122 


2000-1-17 (4711) 


ACIS-S3 


8.6 


607 


1967 


2000-12-07 (5038) 


ACIS-S3 


98.8 


9030 


1044 


2001-4-25 (5176) 


ACIS-S3 


17.8 


1800 


2831 


2001-12-12(5407) 


ACIS-S3 


49.4 


6226 


2832 


2002-5-15 (5561) 


ACIS-S3 


44.3 


6427 


3829 


2002-12-31 (5791) 


ACIS-S3 


49.0 


9277 


3830 


2003-7-8 (5980) 


ACIS-S3 


45.3 


9668 


4614 


2004-1-2 (6157) 


ACIS-S3 


46.5 


11856 


4615 


2004-7-22 (6359) 


ACIS-S3 


48.8 


17979 


4640+4641+5362 


2004-8-26~9-5 


ACIS-S +LETG 


289.0 


16557^ 


+5363+6099^ 


(~6400) 








5579+61782 


2005-1-12(6533) 


ACIS-S3 


48.3 


24939 


5580+63452 


2005-7-14 (6716) 


ACIS-S3 


44.1 


27048 



^ Day since SN. 

2 These observations were splitted by multiple sequences which were combined for the analysis. 

3 Photon statistics are from the zeroth-order data. 



N 



1 arcsecond 




(a) 2000-1 -1 7 (471 1 ) (b) 2001 -4-25 (51 76) (c) 2002-5-1 5 (5561 ) 






(d) 2003-7-8 (5980) (e) 2004-7-22 (6359) (f) 2005-7-1 1/1 6 (671 6) 



Figure 1. The 0.3—8 keV band ACIS images of SNR 1987A. These images have been deconvolved with the detector PSF 
and then smoothed for the purpose of display following methods described in literatures (Burrows et al. 2000; Park et al. 
2002). Darker gray-scales are higher intensities. 
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Figure 2. X-ray spectrum of SNR 1987A as taken with 
the ACIS-S3 in 2000 January and 2005 July. The best-fit 
two-shock models are overlaid on each spectrum. 

3. X-RAY IMAGES 



X-ray images of SNR 1987 A from the ACIS observa- 
tions are presented in Figure 1. Figure 1 shows images 
from six epochs which are separated roughly by ^^1 yeiir 
from each other. The continuous brightening of the X-ray 
emission from the SNR is evident. The significant bright- 
ening of the remnant makes the X-ray morphology now a 
complete ring. This X-ray morphology suggests that the 
blast wave shock is now engulfing the entire inner ring 
rather than small clumps of the dense CSM. 



4. X-RAY FLUX 



Figure 2 shows the ACIS spectrum of SNR 1987 A from 
two epochs of 2000 January and 2005 July. The overall 
increase of the X-ray flux in the last 5.5 years is evident. 
Previous works reported that the electron temperature of 
the X-ray emitting plasma, measured by a single shock 
model, was decreasing from kT ~ 3 keV (day ~4600) 
to kT ~ 2.2 keV (day ^6200) (Park et al. 2004;2005a). 
Our latest data indicate an electron temperature of kT ~ 
1 .6 keV as of 2005 July (day ~6700), which confirms the 
overall spectral softening of SNR 1987 A. As the overall 
spectral shape changes with the photon statistics signifi- 
cantly improving, the single-shock model, however, be- 
came unable to adequately fit the recently obtained data. 
A two-shock model is more useful approximation than 
a single-shock model in order to describe the observed 
spectrum by characteristically representing the deceler- 
ated and the fast shock components (e.g.. Park et al. 
2004). We thus fit the X-ray spectrum observed in each 
eleven epoch with a two-component plane-parallel shock 



model (Borkowski et al. 2001). The detailed description 
of the two-shock spectral modeling of SNR 1987A can 
be found elsewhere (Park et al. 2005c). 

In the spectral fits, we fix some elemental abundances 
(relative to solar) at the values appropriate for the inner 
ring and the LMC ISM because contributions of the line 
emission from those species in the observed energy range 
are negligible: i.e.. He (= 2.57), C (= 0.09) (Lundqvist 
& Fransson 1996), Ca (= 0.34), Ar (= 0.54), and Ni (= 
0.62) (Russell & Dopita 1992). For other species of N (= 
0.76), O (= 0.094), Ne (= 0.29), Mg (= 0.24), Si (=0.28), 
S (= 0.45), and Fe (= 0.16), we use the abundances mea- 
sured with the spectral analysis of our deep LETG obser- 
vations (Zhekov et al. 2005b), because we believe that 
the high resolution dispersed spectrum obtained from the 
deep LETG observations provides the best information 
on the elementral abundances. The best-fit models then 
indicate the electron temperatures of kT = 0.22—0.31 
keV and 2.2—3.2 keV for the soft and the hard compo- 
nents, respectively. The soft component is most likely 
in coUisional ionization equilibrium (ionization timescale 
Hgt 10^^ cm~^ s), while the hard component is in non- 
equilibrium ionization condition {ngt ~ 2 x 10^^ cm~^ 
s). The best-fit total foreground colunon is Nh = 2.35 x 
10^1 cm-2. 

Based on the two-shock model, we present the X-ray light 
curves of SNR 1987A in Figure 3. The soft X-ray light 
curve shows that the observed flux significantly deviates 
from the extrapolation of the simple model used in previ- 
ous works (the short-dashed curve) after day '--^6200. We 
interpret that this up-turn of the soft X-ray flux is caused 
by the shock recently beginning to interact with the en- 
tire inner ring rather than only with small protrusions of 
the CSM. In fact, the modified model (the sohd curve), 
which considers the X-ray flux separately before and af- 
ter the shock interacts with the dense CSM in order for a 
better description of such a condition, can successfuUy fit 
the overaU soft X-ray light curve (Park et al. 2005b). 

The light curves of the fractional contributions from the 
soft and the hard components to the observed 0.5—2 keV 
flux are presented in Figure 4. The contribution from 
the soft component was small and then continues to in- 
crease for the last five years to become dominant since 
day ~6200. This long-term change in the fractional light 
curves suggests that the decelerated portion of the shock 
front by the interaction with the dense CSM has contin- 
uously increased and now the shock front is for the most 
part decelerated by the dense inner ring. The soft X-ray 
light curve and the fractional flux variations thus consis- 
tently support that the blast wave shock began to interact 
with the entire inner ring since day ~6200. 



5. RADLVL EXPANSION RATE 



The radial expansion of the overall X-ray remnant of SN 

1987A has been reported (Park et al. 2002;2004). The 
expansion rate was estimated to be ~4200 km s~^ un- 
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Figure 3. Composite light curves of SNR 1987 A (as presented in Park et al. 2005b). The ROSAT data are taken from 
Hasinger et al. (1996). The radio fluxes were obtainedwith the Australian Telescope Compact Array (ATCA) (provided by 
L. Staveley-Smith) and arbitrarily scaled for the purpose of display. The solid curve is the best-fit model to fit entire soft 
X-ray (0.5—2 keV) light curve (Park et al. 2005b). The short-dashed curve is the best-fit model from Park et al. (2004), 
which is extrapolated after day 6000. The long-dashed line is a linear fit to the ROSAT data (Burrows et al. 2000), which 
is extrapolated to day 7000. The inset is the same light curve plot in a log-scale. 



tU day ~5800 (Park et al. 2004). With the latest data, 
we perform X-ray measurements of the radial expansion 
rate of SNR 1987A using a more sophisticated image 
modeling than that by Park et al. (2002;2004). The de- 
tailed description of our image analysis for the expan- 
sion measurements can be found elsewhere (Racusin et 
al. 2005). We find a "break" in the radial expansion 
rate at around day 6200 where the rate becomes signif- 
icantly lower from ^3800 km to ^1600 km s^^ (Fig- 
ure 5). The deceleration of the radial expansion rate of 
SNR 1987A is perhaps expected considering our picture 
of the blast wave evemtuaUy entering the dense CSM all 
around the inner ring. It is remarkable that the predicted 
deceleration of the expansion rate occurs at day ^^6200 
which is coincident with the results from the soft X-ray 
light curves. The radial expansion rate of SNR 1987A is 
thus self-consistently supportive of the shock beginning 
to interact with the entire inner ring at day ~6200. 



6. DISPERSED SPECTRUM 



The dispersed X-ray spectrum of SNR 1987 A obtained 
with deep LETG observations is presented in Figure 6. 
The fine structures of X-ray hues are clearly resolved 
with good photon statistics. The detailed discussion of 
the results from these LETG data can be found in Zhekov 
et al. (2005a;2005b). We note that we chose the roll- 
angle of the observations in order to align the dispersion 
axis in the north-south direction. With this observation 
setup, assuming that the bulk of X-ray emission origi- 
nates from the interaction between the blast wave and the 
inner ring, we expected that the Doppler shifts due to the 
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Figure 4. The 0.5—2 keV band fractional light curves of 
SNR 1987 A based on the two-shock model (as presented 
inParketal. 2005b). 
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Figure 5. The radial expansion rate of SNR 1987 A. 
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Figure 6. The first-order ( +1, solid line) and the zeroth- 
order (dashed line) spectrum of SNR 1987 A from the 
LEGT observations ( as presented in Zhekov et al. 2005a ). 
Some identified lines from key elemental species are 
marked. 

the orientation of the inner ring (an inchnation angle of 
~45° with the north side toward the Sun) would result 
in systematic distortions of the dispersed images of SNR 
1987A: i.e., the overall SNR image should be compressed 
in the negative (to = — 1) arm while it is stretched in the 
positive (to = +1) arm. The measured line widths indeed 
reveal systematically larger broadenings in the positive 
arm than those measured in the negative arm (Figure 7). 
These line profiles provide the first direct evidence that 
the X-ray emission is originating from a "plane" contain- 
ing the inner ring rather than from a spherical volume. 
These line broadenings imply shock velocities of v ~ 
340-1700 km s-^ 

The deep LETG observations allow us to directly mea- 
sure the line flux ratios from various X-ray lines of SNR 
1987A (Zhekov et al. 2005a). In Figure 8, we present flie 
ranges of the electron temperature (kT) and the ioniza- 
tion timescale {riiA) derived from the plane -parallel shock 
model, which agree with the measured "G-ratios" of the 
Hea triplets (G = [f + i]/r, where f, i, and r are the 
forbidden, intercombination, and resonance line intensi- 
ties, respectively) and the Hea/Lya ratios from O and Si. 
Broad ranges of the electron temperature {kT ~ 0.1—2 
keV) and the ionization timescale {Uet ^ 10^^ — 10^^ 
cm~^ s) are consistent with the measured Une ratios. 



7. DISCUSSION 
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Figure 7. Measured line widths (FWHM) for the 
(diamonds) and the "—1 " (squares) LETG arms (as pre- 
sented in Zhekov et al. 2005a). The solid curve presents 
the resolving power of the LETG. The dashed and dot- 
ted curves are best-fit line-broadening parameters for the 
cases with and without shock stratification, respectively 
(see Zhekov etal. 2005a for the details). 

ported by several observations. 

(1) A simple model assuming a single shock propagating 
into ambient medium with an exponential density profile 
has successfully described the soft X-ray light curve for 
the previous ~13 years. However, an up-turn of the soft 
X-ray flux since day ~6200 cannot be described with the 
simple constant- velocity shock model. A significant con- 
tribution from the decelerated shock needs to be consid- 
ered in order to adequately fit the data (Park et al. 2005b). 

(2) The fractional contribution to the observed 0.5—2 ke V 
flux from the soft component of the two-shock model 
has continuously increased. The soft component repre- 
sents the X-ray emission from the decelerated shock and 
becomes dominant since days ^6000—6200 (Park et al. 
2005b). (3) The X-ray radial expansion rate significantly 
reduces since day ~6200 (Racusin et al. 2005). (4) 
The shock kinematics obtained by the LETG observa- 
tions (at day ^6400) indicate velocities of ^340—1700 
km s~^ for the X-ray emitting plasma (Zhekov et al. 
2005a). These velocities are significantly lower than 
those estimated from earlier X-ray and radio images (v ~ 
3000-4000 km s'^). (5) The soft X-ray images indicate 
that the brightening has become global since day ^6000 
(Park et al. 2005b). (6) The optically bright spots become 
prevailing over the entire iimer ring by day ~6000 (e.g., 
Plate 1 in McCray 2005). 



The bulk of soft X-ray emission of SNR 1987A most 
likely originates from the interaction between the blast 
wave shock and the clumpy structure of the pre-existing 
(before the SN) dense equatorial CSM. The emergence 
of the first optical spot in day ^-^3700 indicated the be- 
ginning of such interactions. Days ~6000— 6200 appear 
to mark another important milestone in the evolution of 
SNR 1987A: i.e., the blast wave reaches the dense CSM 
around the entire inner ring. This interpretation is sup- 



On the other hand, the hard X-ray light curve shows sig- 
nificantly lower flux increase rate (Figure 3). It is remark- 
able that the hard X-ray light curve is similar to that of the 
radio data. The radio emission appears to primarily orig- 
inate from the synchrotron emission behind the reverse 
shock (Manchester et al. 2005). The similarity between 
the radio and hard X-ray light curves might suggest that 
the hard X-rays originate from the same population of 
electrons to produce radio synchrotron emission. How- 
ever, a better correlation of the radio image with the hard 
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Figure 8. Allowed kT vs. Ugt ranges to match the measured G-ratios (dotted curves) and Hea/Lya ratios (dashed curves) 
of O (left panel) and Si (right panel) lines from SNR 1987A (as presented in Zhekov et al. 2005a). The three curves for 
each ratio correspond to the measurement and ±la uncertainties, respectively. 



X-ray image than with the soft X-ray image is not clear 
(Park et al. 2005b). Thus, the lower flux increase rate in 
the hard X-ray band might simply be caused by the con- 
tinuous softening of the X-ray spectrum due to the shock 
interaction with the dense CSM. Follow-up monitoring of 
the light curve, the morphology, and the spectral proper- 
ties in the hard X-ray band will be essential to unveil the 
true origin of the hard X-ray emission of SNR 1987A. 

With the recent deep LETG observations, we for the first 
time measure the individual X-ray line profiles from SNR 
1987A (Zhekov et al. 2005a). The systematic difference 
in the line broadening between the positive and negative 
dispersion arms provide direct observational evidence for 
the X-ray emission originating in the equatorial plane of 
the inner ring. More surprisingly, the derived velocity 
(v ~ 340—1700 km s~^) of the X-ray emitting plasma 
is significantly lower than the average shock velocities 
estimated with the previous X-ray and radio data (v ^ 
3000—4000 km s^^). The low velocity suggests that the 
blast wave has been significantly decelerated by the dense 
inner ring, which is consistent with the results from the 
X-ray light curve and image analyses. 

The X-ray line ratio measurements indicate that the X- 
ray emitting plasma of SNR 1987 A cannot be adequately 
described by a single combination of the electron temper- 
ature and the ionization timescale. Instead, the observed 
X-ray lines originate from wide ranges of the plasma tem- 
perature and the ionization age. The multi-phases of the 
X-ray emitting thermal plasma of SNR 1987A was sug- 
gested by the undispersed spectrum and is clearly con- 
firmed by the dispersed spectrum. This plasma structure 
is most likely due to the interaction of the blast wave with 
the complex density gradients near the boundary of the 
dense inner ring. 
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